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Abstract 
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. We investigate K and K condensation in /3-equilibrated hyperonic matter 

within a density dependent hadron field theoretical model. In this model, 
\ baryon-baryon and (anti)kaon-baryon interactions are mediated by the ex- 

change of mesons. Density dependent meson-baryon coupling constants are 
obtained from microscopic Dirac Brueckner calculations using Groningen and 
Bonn A nucleon-nucleon potential. It is found that the threshold of antikaon 
condensation is not only sensitive to the equation of state but also to antikaon 
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CN \ optical potential depth. Only for large values of antikaon optical potential 



in 



depth, K condensation sets in even in the presence of negatively charged 



I/"") ■ hyperons. The threshold of K° condensation is always reached after K~ con- 

densation. Antikaon condensation makes the equation of state softer thus 
, resulting in smaller maximum mass stars compared with the case without 

^ (— I any condensate. 

Oh| PACS: 26.60.+C, 21.65.+f, 97.60.Jd, 95.30.Cq 
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I. INTRODUCTION 

a 

The composition and structure of neutron stars depend on the nature of strong interac- 
tion. Neutron star matter encompasses a wide range of densities, from the density of iron 
nucleus at the surface of the star to several times normal nuclear matter density in the core. 
Since the chemical potentials of nucleons and leptons increase rapidly with density in the 
interior of neutron stars, several novel phases of matter with large strangeness fraction such 
as, hyperonic matter, condensates of strange mesons and quark matter may appear there 

I- " ' 

It was first demonstrated by Kaplan and Nelson within a chiral SU(3)l X SU(3)r model 
that K~ meson may undergo Bose-Einstein condensation in dense matter formed in heavy 
ion collisions 0. In this model baryons directly couple with (anti)kaons. The effective mass 
of antikaons decreases with increasing density because of the strongly attractive i^~-baryon 
interaction in dense matter. Consequently, the in-medium energy of K~ mesons in the zero- 
momentum state also decreases with density. The s-wave K~ condensation sets in when the 
energy of K~ mesons equals to its chemical potential. Later, K~ condensation in the core 
of neutron stars was studied by other groups using chiral models JJI. 

Also, Bose-Einstein condensation of K~ mesons was investigated in the traditional me- 
son exchange picture known as relativistic mean field (RMF) model [|],[5]||. Within the 
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framework of RMF model, baryon-baryon and (anti)kaon-baryon interactions are treated in 
the same footing i.e. those are mediated by the exchange of mesons It was noted in 

all these calculations that the typical threshold density of K~ condensation in nucleons-only 
neutron star matter was about 2-4no, where n is normal nuclear matter density. However, 
the threshold of antikaon condensation is sensitive to the antikaon optical potential and 
depends more strongly on the equation of state. With further inclusion of hyperons, K~ 
condensation was found to occur at higher densities ||,[5],|6],|8]|| . Recently, we have studied 
K° condensation along with K~ condensation in neutron stars using a relativistic mean 
field model and its influence on the structure of neutron stars The threshold density 

of K° condensation is always higher than that of K~ condensation. Antikaon condensate 
makes the equation of state (EoS) softer, thus resulting in smaller maximum mass stars 
compared with the case without any condensate. Employing the EoS including both K~ 
and K° condensates, it was predicted that a stable sequence of superdense stars called 
the third family branch [TU] might exist beyond the neutron star branch. The compact stars 
in the third family branch have smaller radii than those of the neutron star branch [[|. 

Besides EoS and antikaon optical potential depth, the threshold density of antikaon con- 
densation is very much sensitive to the behaviour of antikaon energy and electron chemical 
potential at high density. The role of nucleon-nucleon and (anti)kaon-nucleon correlation on 
antikaon condensation were investigated by Pandharipande and collaborators [0,12^. They 
found that strong nucleon-nucleon and (anti)kaon-nucleon correlation raised the critical 
density for antikaon condensation to higher densities and predicted antikaon condensation 
might not be a possibility in neutron stars | 12fl . The electron chemical potential used in the 
above mentioned calculations was obtained from modern realistic nucleon-nucleon interac- 



tions |L3],[14 



In this work, we are interested to find out how many body correlations which may be 
taken into account by density dependent meson-baryon couplings in a relativistic field theo- 
retical model, affect the threshold of antikoan condensation in neutron star matter. There is 
a growing interest to derive a quantum hadron field theory from a microscopic approach to 
nuclear interactions. The motivation of such an approach is not only to retain the essential 
features of quantum hadrodynamics (QHD) but also to deal with the complicated many 
body dynamics of strong interactions |T6|JT7| , p!8| . An appropriate and successful microscopic 
approach to in medium nuclear interactions follows from Dirac Brueckner (DB) calcula- 
tions. Various groups performed DB calculations with realistic nucleon-nucleon interactions 
and reproduced empirical saturation properties of symmetric nuclear matter reasonably well 
|I^ , |20|pT| , |22|p3| , |2^p5| , |26[| . Also nucleons-only neutron star matter was calculated in Dirac 
Brueckner theory [^,^,^]. It is worth mentioning here that the bulk of the screening of 
nucleon-nucleon interaction in medium is taken into account by the local baryon density 
dependent DB self energies ||17|| . This makes relativistic many body dynamics to be approx- 
imated by a density dependent relativistic hadron (DDRH) field theory fl7| , [T8|j . A covariant 
and thermodynamically consistent DDRH field theory is obtained by making interaction 
vertices as Lorentz scalar functionals of baryon field operators. In the mean field approx- 
imation this model reduces to the relativistic Hartree description with density dependent 
meson-nucleon couplings. The density dependent meson-nucleon couplings are obtained 
from Dirac Brueckner self energies calculated with Bonn, Groningen and phenomenologi- 
cal density dependent potentials [^,[JT],[^,[£]|,[£|,[^] . The variational derivatives of vertices 
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with respect to baryon fields give rise to rearrangement terms in baryon field equations |18| 



Brockmann and Toki [T(| first applied DDRH model without rearrangement terms to study 
finite nuclei. Recently, DDRH model with rearrangement terms has been exploited to inves- 
tigate deformed nuclei pi] , hypernuclei [33], asymmetric nuclear matter and exotic nuclei 
34] , |3"5|| and neutron star properties |36| . 

In this paper, we investigate antikaon condensation in beta equilibrated hyperon matter 
relevant to neutron stars and its role on the composition and structure of the compact stars 
in DDRH model. The paper is structured in the following way. In Sec. II, we describe 
DDRH model and different phases of matter. Parameters of the model are discussed in Sec. 
III. Results of our calculation are explained in Sec. IV. And Sec. V provides a summary 
and conclusions. 



II. FORMALISM 

Here, we discuss a phase transition from hadronic matter to antikaon condensed phase 
in compact stars. The hadronic phase is described within the framework of DDRH model. 
This phase is composed of all species of the baryon octet, electrons and muons. Therefore, 
the total Lagrangian density in the hadronic phase is written as C = Lb + A- m DDRH 
model, baryon-baryon interaction is given by the Lagrangian density (£b) p6| , 



C B = J2^b [ Hud* 1 -m B + 9aBO - QuBl^ - -g P BlfjTB • P M + ^sbTb • S ) ^ B 



B 



and 



+ l - (dpaVa - mla 2 ) + l - (d^dPS - m 2 5 2 ) - 

+~m^y^ - -p^ ■ p^ + \m 2 pPll ■ (I) 



A = E^(*7^-m,)Vi • (2) 
l 

The field strength tensors for vector mesons are given by 

fT = d»p v - d u ff (3) 

Here \& b denotes the isospin multiplets for baryons and the sum goes over baryon multi- 
plets B = N,A,E,H; ipi {I = e , A 4 ) is lepton spinor and tb is isospin operator. The interactions 
among baryons are mediated by the exchange of a, uj and p mesons. In addition to these 
mesons the scalar-isovector meson S is also included. And this is important for an asymmet- 
ric system. Though the structure of DDRH Lagrangian density closely follows that of RMF 
model, there are important differences between those models. In RMF calculations with 
density independent meson-baryon coupling constants, non-linear self interaction terms for 
scalar and vector fields are inserted to account for higher order density dependent contri- 
butions. But this is not needed here as meson-baryon vertices g a B, where a denotes a, u>, 
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p and 5 fields, are dependent on Lorentz scalar functionals of baryon field operators and 
adjusted to the Dirac-Brueckner-Hartree-Fock (DBHF) calculations [34.361 . 



There are two choices for the density dependence of meson-baryon couplings. One is the 
scalar density dependence (SDD) and the other one is the vector density dependence (VDD) 



34[. Here we consider meson-baryon couplings g a B(f>) to depend on vector density because 



it gives a more natural connection to the parameterization of DB vertices. For VDD case, 



the density operator p has the form, p=y where j M = ^-y^. 

As vertices g a BS are Lorentz scalar functional of baryon field operators, the variation of 
C with respect to ^I>b gives, 

+ i^r^r ( 4 ) 



o4 B dQ B 9 fa S^b 

The rearrangement term E^ r ) = J2b J^H^" which originates from the second term of Eq. 
(4), naturally introduces additional contribution to vector self-energy |l7| , p!8| , p4| , p6| . This is 
an important difference between RMF and DDRH theory. 

Here, we perform our calculation in the mean field approximation (MFA). In this ap- 



proximation vertex functionals are reduced to simpler forms using Wick's theorem |34| , |37 
The operator p is replaced by ground state expectation value p i.e. < p >= p. Hence 
meson-baryon vertices become function of total baryon density in the hadronic phase, 

< 9aB(p) >= 9aB(< P >) = 9ob(p)- (5) 

This is known as vector density dependence of vertices |T7|J54] , |3^[] . In MFA adopted here, 
meson fields are replaced by their expectation values. Only the time-like components of 
vector fields, and isospin 3-components of p and 5 fields have non-vanishing values in a 
uniform and static matter. The mean meson fields are denoted by a, uj , p 03 and 5. Therefore, 
the meson field equations in the hadronic (h) phase are given by 

(6) 
(7) 

(8) 

(9) 

where r 3 b is the isospin projection of baryon b = n, p, A, E~, S°, S + , E~, S° and scalar and 
vector densities of baryon b in the hadronic phase are 



m 2 a a 


\ -« h.s 

= }^9abn b ' , 

6 




= J2 9^n h h , 
b 


m 2 p po3 


= \ Y.9 P br-ibn 1 l 
z b 


m 2 S 


= 9 S 95bT3 b n h b ' 
L b 



h,s / 7 / ^ 2J 6 + 1 fkF„ ml 



ml 



■ kF h + \ k F , + mt 2 
[k Fb Jk Fb 2 + mf - mfln V * b -\ , (10) 



n h b = <^ b ^ b >=^ . (11) 
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Here Jb is the spin projection of baryon b. The rearrangement self energy modifies the 
baryon field equation compared to the RMF case [15| . 



[ 7 „ - E°,) - (m 6 - Efo)]V* = 0. (12) 

Here ipj, is the Dirac spinor for baryon b. The total vector self energy for baryon b in the 
hadronic phase is 

v o _ v°(°) _l v°( r ) n r>\ 

^b,h — ^b,h ' ■ i 16 ) 

Now the usual vector self energy is 



S bS? = 9ubU0 + -g p bT3bP03 ■ (14) 



Also, the rearrangement term which is the second term in Eq. (4), simplifies to |34 



v ow v^r dg a b h,s . 9g^b h , Idgpb h ^-dgsb x h S i / 1f -\ 

= D"^™* + + tQ^TnfW ~ 2^ 6n » ] > ( 15 ) 

Similarly, the expression of scalar self energy for baryon b is given by 

^t,h = g*bv + ^gsbTsbS. (16) 

One can immediately define effective baryon mass as ml = mb — ^lh- And this differs for 
members of isospin multiplets due to S meson. 

To obtain the EoS (pressure versus energy density) of the pure hadronic phase, the 
equations of motion for mesons and baryons (Eqs. (6)-(9) and Eq. (12)) are solved self- 
consistently along with effective baryon mass in the mean field approximation keeping into 
consideration other constraints such as baryon number conservation, charge neutrality and 
beta equilibrium. The system is charge neutral and the condition of /3-equilibrium is main- 
tained. The charge neutrality condition is 

Q h = Y,to n b -n e -7v = 0, (17) 

b 

where qb and are electric charge and the number density of baryon b in the pure hadronic 
phase, respectively and n e and are number densities of electrons and muons respectively. 
In compact star interior, chemical equilibrium is maintained through weak interactions such 
as B\ — > B 2 + l + i>i and B 2 + 1 — > B\ + ui where B\ and B 2 are baryons and 1 stands for 
leptons. Therefore the generic equation relating chemical potentials for the above mentioned 
generalised /3-decay processes is 

Hi = bi/j, n - qt^e , (18) 

Here 6« and g» are the baryon number and charge of ith baryon and fi n and /i e are the 
chemical potentials of neutron and electron, respectively. The chemical potential of baryon 
b in the hadronic phase is expressed as 
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Ai6 = V*Si + ™J2 + E$>+E? r >. (19) 

It is noted that unlike RMF model, the rearrangement term appears in the expression of 
baryon chemical potential in DDRH model. In neutron stars, electrons are converted to 
muons by e~ — > /i~ + + z/ e when the electron chemical potential becomes equal to the 
muon mass. Therefore, we have fi e = /x M in compact stars. Equation (18) implies that there 
are two independent chemical potentials \i n and \i e corresponding to two conserved charges 
i.e. baryon number and electric charge. The energy density (e h ) is related to the pressure 
(P h ) in this phase through the Gibbs-Duhem relation 

P h = Y,Vi n i- £h ■ (2°) 

Here fa and rii are chemical potential and number density for i-th species. The expression 



of energy density in the hadronic phase is [36 



1111 

h 2 2 , L 2 2 i 1 2 2 , L 2x2 

£ = 2 m <x a + 2 m - w o + 2 m P Po3 + 2 m * S 

2J h + i rH,.-, *om/o.q ^ l f Kp i 



+ £ I ^ ( k2 + < 2 ) 1/2 k 2 dk + Y.-J V + W) V2 ^ 2 dfe- (21) 

k 27T Z JO j 7T Z JO 

The rearrangement term does not contribute to the energy density explicitly, whereas it 
occurs in the pressure through baryon chemical potential. It is the rearrangement term 
that accounts for the energy-momentum conservation and thermodynamic consistency of 



the system fl8fl 



The pure antikaon (K) condensed phase is composed of baryons, leptons and antikaons 
which are in chemical equilibrium under weak interactions and maintain local charge neu- 
trality. The baryon-baryon interactions here are described by the Lagrangian density of 
DDRH model. It is worth mentioning here that the meson-baryon couplings depend on the 
total baryon density in this phase. In this phase, baryons are embedded in antikaon conden- 
sates. Earlier it was noted that baryons in the pure hadronic and antikaon condensed phase 
behaved differently because of their dynamical nature [0,H|§. It was attributed to different 
mean fields which baryons experienced in those pure phases. We adopt here a relativistic 
field theoretical approach for the description of (anti)kaon-baryon interaction In this 

model (anti)kaon-baryon interactions are mediated by a, u>, p and S meson. The Lagrangian 
density for (anti)kaon interaction in the minimal coupling scheme is, 

C K = D*KD»K - m§KK , (22) 

where the covariant derivative = <9 M + ig w K^^ + +W p k t k • P M /2. The isospin doublet 
for kaons is denoted by K = (K + ,K°) and that for antikaons is K = (K~,K°). It is to 
be noted that the coupling constants of (anti)kaon-baryon interactions are considered to be 
density independent. The effective mass of (anti)kaons in this minimal coupling scheme is 
given by 

m* K = m K - g uK a - ^-gsxr^ , (23) 
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where mx is the bare kaon mass. Here also the effective mass of K~ and K° differ due 
to the inclusion of the scalar-isovector 5 meson. The dispersion relation representing the 
in-medium energies of K = (K~ , K°) for s-wave (k = 0) condensation is given by 

VK- t K° = m *K- 9wKU =F ~g p KP03 , (24) 

where the isospin projection r 3 # = =Fl for the mesons K~ (— sign) and K° (+ sign) are 
explicitly written in the expression. For s-wave condensation, densities of antikaons are 
given by 

n K-,K° = 2 (^k-,k° + 9uKW ± -g P KPoz \ KK = 2m* K KK . (25) 

In the mean field approximation, the meson field equations in the presence of antikaon 
condensates are given by 



ml<r = 9abnf ,s + g uK ^2 n K > (26) 

K 



m l u o = 9^ - 9ojK Y riK , (27) 

b R 
1 - 1 

m lP03 = T £ 9 P bT3 b nb + -9pK T 3K n K , (28) 
b K 

m s s = gsbTsbn*' 8 + \g& K T 3i<nK , (29) 
1 b 1 R 

where nf ,s and n§ are scalar and vector density of baryon b in the antikaon condensed 
phase and have the same forms as in Eqs. (10) and (11). The meson field equations here 
remain same in structure as RMF ones [|7]||ll, ^° u ^ ^ ne constant meson-baryon couplings are 
replaced by their density dependent counterparts. 

The total energy density and pressure in the antikaon condensed phase are given by ||J| 

£ K = \ m l° 2 + \ m l u l + \ m lPm + \ m ^ 2 

+ Y 2 A ± T- [ kFb (k" + mfy /2 k 2 dk + Y-2 [ KF '(k 2 + mWk 2 dk 
^ 2ir z Jo j ir z Jo 

+m* K (n K - + n^o) , (30) 

and 

r>K 1 22,1 22,1 2 2 ^ 2jr2 , v^Q(r) K 

p = ~2 m ^ + 2 " o + 2 mpPo3 ~ 2 m& * 2^ n b 

1 ^ 2J b + 1 rk Fb fc 4 dk I v — [ KFl fc4 dk 

+ 3 V 2vr 2 Jo {k 2 + mfy/ 2 + 3^^ 2 Jo {k 2 + m 2 ) 1 / 2 ' l>! 

where ^% is the rearrangement term in the antikaon condensed phase and has the same 
form as in Eq. (15), but all quantities in the equation are replaced by the corresponding 
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quantities of the antikaon condensed phase. Since antikaons form s-wave condensates, they 
do not contribute to the pressure directly. Actually the effect of antikaons in the pressure 
term comes through the meson fields. 

In the core of neutron stars, various strangeness changing processes such as iV ^ iV + K 
and e~ # K~ + z/ e may occur Here N = (n,p) and K = (K~, K°) denote the isospin 

doublets for nucleons and antikaons, respectively. From the above reactions in chemical 
equilibrium, we obtain the conditions for antikaon condensation 

— = /Aff- = i (32) 

/^o = , (33) 

where fix- and fi^o are respectively the chemical potentials of K~ and K°. The charge 
neutrality condition in the antikaon condensed phase is 

Q K = ^ n b ~ n K- - n e - = 0. (34) 

b 

It was noted in RMF model calculations that antikaon condensation could be either first 
order or second order phase transition depending on the parameter set of the model and 
antikaon optical potential depth If the phase transition is of first order, the mixed 

phase is to be determined by the Gibbs conditions and global baryon and electric charge 
conversation laws because we have conserved baryon and electric charges represented by two 
chemical potentials \i n and /i e [BHJ. The Gibbs phase rules read, 



P h = P R , (35) 
A*» = A*f , (36) 

where /z£ and are chemical potentials of baryon b in the pure hadronic and K~ con- 
densed phase, respectively. The conditions of global charge neutrality and baryon number 
conservation are imposed through the relations 

(i - x)Q h + xQ R = 0, (37) 

n b = (1 - xK + Xnf , (38) 

where x is the volume fraction of K~ condensed phase in the mixed phase. The total energy 
density in the mixed phase is 

e=(l- x )e h + X e R . (39) 



III. PARAMETERS 

A. Meson-nucleon couplings 

In DDRH model, the dependence of meson-nucleon vertices on total baryon density is 
obtained from microscopic DB calculations of symmetric and asymmetric nuclear matter. 
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The density dependent vertices in the RMF model are related to DB self energies in the local 
density approximation • Equating self-energies of infinite nuclear matter in RMF and 



DB calculations, we obtain 



^RMF i s^DB 



= E , (40) 

where <p a represents the field for a meson. Putting the value of <p a as given by the meson 
field equations in the presence of nucleons, the above relation simplifies to 



2 Z DB = g 2 a n a (41) 



where n a is the density corresponding to 4> a field. Scalar and vector self energies for neutrons 
and protons were obtained in DB calculations of asymmetric nuclear matter using Gronin- 
gen potential [p6| , p9| , ^0| . Using these values of E DB , one immediately obtains the density 
dependent meson-baryon couplings [[34J. A suitable parameterization for density dependent 
couplings was made in Ref. [j34]1 . It has the form 



1 + b a (p p + d a ) 
9M = a*— — 1 —, — - — 42 

1 + c a {p/po + e a ) 2 

where po — 0.16/m~ 3 and parameters of the fit are listed in Table I of Ref. |j34 |. However, 
the results of infinite nuclear matter calculation in DDRH model using the above mentioned 
parameterization deviated from those of DB calculations |34| because momentum dependent 



DB self energies were mapped onto the momentum independent DDRH self energies. There- 
fore, momentum dependent vertices with the additional constraint that the energy density 
in DB and DDRH are same i.e. e DB = £ DDRH was proposed ||34|| . Momentum corrected 
meson-nucleon vertices are given by 



g a {k F ) = g a (k F )\jl + ( a k 2 F = ( a {k F )g a (k F ). (43) 

Momentum corrections Cr = 0- 00804 fm 2 and (^=0. 00103 fm 2 to u-nucleon and tu-nucleon 
were obtained from DB calculations of symmetric nuclear matter Using the momen- 



tum corrected vertices, Dirac Brueckner EoS for symmetric matter was reproduced well [34 



Similarly, momentum correction ( p to p-nucleon vertex was calculated from DB calculations 



of neutron matter p4[ . This correction was inserted in the DB self energies and the mo- 
mentum corrected p-nucleon vertex was calculated. Later the density dependence of the 
momentum corrected p-nucleon vertex was parameterized using Eq. (42). The parameters 
of this fit are given by Table II in Ref. |34j| . We adopt this parameterization of density 



dependent couplings and the momentum correction prescription in our calculation. Also, 
we denote this as Groningen parameter set. In Table I, we show meson-baryon couplings for 
Groningen set at saturation density (n = 0.18/mr 3 ). The momentum correction modifies 
the rearrangement term as is to be replaced by this is given by f34| 



dg aB {k F ) ~ n ,dg aB (k F ) ( a k F 2 g aB (k F ) 

__ = aM __ + __. (44 ) 

In this calculation, we also exploit density dependent meson-nucleon vertices obtained 
from DB calculations using Bonn A potential. The parameterization of vertices is taken from 
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Ref. |30]]. This parameter set is denoted as Bonn A parameter set. For Bonn A set, p meson- 
nucleon coupling is chosen as a constant. Also, 5 meson is not taken into consideration for 
Bonn A potential. Meson-nucleon coupling constants at saturation density no = 0.159/m -3 
are listed in Table I. 



B. Meson-hyperon couplings 

In the absence of DB calculation including hyperons, density dependence of meson- 
hyperon vertices are obtained from density dependent meson-nucleon couplings using hy- 
pernuclei data || and scaling law |j33| . This scaling law states that the self energies and 



vertices of hyperons and nucleons are related to each other by their free space coupling 



constants g a y and g a N p3|j3q| . 



R aY = = |^ « (45) 

9aN ^aN 9aN 

In RMF model, vector meson-hyperon coupling constants were determined from scaling fac- 



tors obtained from SU(6) symmetry relations of the quark model [§,0. Another possibility 
is to exploit scaling factors calculated in microscopic calculations. However, there is only 
one microscopically derived free space scaling factor R„a =0.49 in the literature [|42[j . We 
use this value in our calculation. Also, we obtain the scaling factors for vector and isovector 
mesons from SU(6) symmetry relations ||, 



1 1 

— 9uE — -^QujNi 

^9 P s = 9 P s = g P N\ 9 p a = 0, 



^flte = gss = 9sn; 98k = 0. (46) 

We obtain g^A and scalar meson couplings to other hyperons from the potential depths 
of hyperons in normal nuclear matter. The hyperon potentials in saturated nuclear matter 
are obtained from the experimental data for the single particle spectra of hypernuclei. In 
DDRH model, the potential depth of a hyperon (Y) in saturated nuclear matter is given by 

UP = Ef J + E% r) - Ey , (47) 

where Ey -* = g^y^o, E y = g^a and E 1 ^ is the rearrangement contribution of nucleons. 
In this calculation, the value of A potential in normal nuclear matter is taken as -30 MeV 



41| , |43[| and that of H is -18 MeV |44|J45|] . The most updated analysis of E atomic data f46 



and other experimental data [|7|| predict a repulsive E-nucleus potential depth. Therefore, 



we adopt a E well depth of 30 MeV in this calculation j|7|]. We find that E hyperons are 



excluded from the system because of this repulsive potential. The scaling factors of A and 
H for Groningen and Bonn A potential are listed in Table II. 

From Table II, we observe that -R^a is 0.4911 corresponding to R a \ = 0.49 for Bonn 
A set. In this case, -R^a is obtained from the A potential depth (U^) as discussed above. 
On the other hand, Keil et al. [R3] obtained a value of R^a = 0.553 for the same value 



of R a A from the x 2 distribution for the deviation of DDRH A single particle energies and 
hypernuclear data. We perform calculations for both the values of R^a for Bonn A set. 
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C. Meson-(anti)kaon coupling constants 



Finally, we need to determine the parameter set for meson-(anti)kaon interactions. Here, 
we do not attribute any density dependence to the vertices of meson-(anti)kaons. The vector 
coupling constants are derived from quark model and isospin counting rule so that 

9uK = 7^9 uN and g pK = 9 p n ■ (48) 

The values of meson-nucleon coupling constants are taken at normal nuclear matter density 
and those are given by Table I. The scalar coupling constant is obtained from the real part 
of K~ optical potential at normal nuclear matter density 

Ur (n ) = -gwKUJo - 9aK(? + S^ r) . (49) 

The scalar isovector 5 meson also couples with (anti)kaons. The coupling of 5 meson with 
(anti)kaons is obtained from simple quark model and this is given by gsx = gsN- The value 
of gsN is obtained from Table I. 

There are experimental evidences that antikaons experience an attractive interaction 
whereas kaons feel a repulsive interaction in nuclear matter 0,^]. It is the depth of 



antikaon optical potential which is an important input in our calculation. The real part 
of antikaon optical potential at normal nuclear matter density was evaluated in coupled 
channel model |50|J51[| and self consistent calculations |52|J53| , |5^] . These model calculations 
give a wide range of values from -120 MeV to -40 MeV for for Ur at n Q . Recently, a 
combined chiral analysis of K~ atomic and K~p scattering data lead to a shallow attractive 



Uk(tio) of -55 MeV |]55| . On the other hand, the analysis of K~ atomic data in the hybrid 
model yielded Ug(no) = —180 ± 20MeV. Therefore, there is no consensus among the 
phenomenological and microscopic potentials both in terms of depth and x 2 values from the 
fits to kaonic atom data. The coupling constants for kaons with cr-meson, g a K, for a set 
values of Ug from -120 MeV to -180 MeV at saturation density for Groningen and Bonn A 
potential are listed in Table III. The cx-K coupling constants for Bonn A set are found to be 
larger than that of of Groningen set. This stems mainly from the smaller g w x^o value for 
Bonn A set compared with that of Groningen set. 



IV. RESULTS AND DISCUSSION 

Here we report the results of our calculation in DDRH model using Groningen set. We 
perform this calculation for antikaon optical potential U^uq) = -120 to -180 MeV. There 
is no K condensation as a first order phase transition for Groningen set and various values 
of Ux{no). Rather, K~ and K° condensation are second order phase transitions in this 
calculation. In this situation, the conditions of antikaon condensation are given by Eqs. 
(32) and (33). Earlier it was found in RMF calculations that antikaon condensation could 
be a second order phase transition depending on the antikoan optical potential and coupling 
constants of the models |7|,|]|| • The threshold densities of K~ and K° condensation in 
/3-equilibrated matter containing n,p,A and leptons for Ur(tiq) = -120 to -180 MeV are 
recorded in Table IV. Besides A hyperons, we also include other species of hyperons into our 
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calculation. However, E hyperons do not appear because of a repulsive S-baryon interaction. 
In Table IV, the critical densities of K condensation in /^-equilibrated n,p,A,S and lepton 
matter are given in the parentheses. The early appearance of hyperons might have important 
effect on the threshold densities of K condensation because hyperons make the equation of 
state soft. It was shown in RMF model calculations that the onset of K condensation was 
delayed to higher densities due to hyperons f|,|^f||8||| . Also, negatively charged hyperons 
diminish the electron chemical potential delaying the onset of K~ condensation. In this 
DDRH model calculation with Groningen set, A hyperons appear first. Consequently, the 
threshold densities of K condensation are shifted to higher densities compared with those 
in nucleons-only matter. With further appearance of negatively charged S~ hyperons, K~ 
condensation occurs at higher densities as it is evident from the values in the parentheses 
in Table IV. For Ux{no) = - 120 MeV, the early appearance of S~ completely blocks the 
onset of both K~ and K° condensation even in the highest density (8n ) considered in this 
calculation. On the other hand, the impact of S~ hyperons on the threshold densities of K° 
condensation for \U^(n )\ > 160 MeV is negligible. This may be attributed to the fact that 
the density of H~ hyperons falls after the onset of K~ condensation. This becomes evident 
when we discuss the particle density graphs in the following paragraphs. From Table IV, 
we note that the threshold density of K condensation shifts towards lower density as the 
strength of |£^(no)| increases. This indicates the threshold of K condensation is not only 
dependent on the EoS, but also sensitive to antikaon optical potential depth. For all values 
of L^j(no), we observe K~ condensation occurs before K° condensation. 

The composition of neutron star matter containing nucleons (n,p), A hyperon, electron 
(e~), muon and K~ and K° mesons for Groningen set and Ux(no) = -160 MeV is 
presented in Fig. 1. Before the onset of K~ condensation, the charge neutrality is maintained 
by protons, electrons and muons in the hadronic phase. We see that A hyperon is the first 
strange baryon to appear in the hadronic phase at 1.99no, where uq = 0.18/m~ 3 . Its density 
rises fast at the cost of neutrons. In this calculation, K~ condensation sets in at 2.65rio- 
As soon as K~ condensate appears, it rapidly grows and readily replaces e~ and \i~ . This 
behaviour is quite expected, as K~ mesons, being bosons, condense in the lowest energy 
state and are therefore energetically favorable to maintain charge neutrality of the system. 
The electron fraction depletes around 4.8n and the proton density becomes equal to that of 
K~ condensate. The appearance of K° condensate is delayed to 6.16no- With the onset of 
K° condensation the abundances of n,p,i^~ and K° become identical leading to an isospin 
saturated symmetric matter This may be attributed to the fact that there is a 

competition in the production of p-K~ and n-K° pairs resulting in a symmetric matter of 
nucleons and antikaons [0,0. Here the system is dominated by A hyperons at high density. 
It is worth mentioning here that the results of DDRH model for Groningen set and Ux{no) = 
-160 MeV resemble those of the RMF model for GM1 set and Ux(rio) = -160 MeV. However, 
antikaon condensation in the latter case was a first order phase transition. 

Besides A hyperon formation, we also consider the role of other hyperons, such as S° and 
H~ on antikaon condensation. In Fig. 2, we note that negatively charged S~ hyperons start 
populating the system at 2.24no, soon after the appearance of A. This further postpones K~ 
condensation to 3.20no- Lepton fractions begin to fall with the onset of negatively charged 
E~. This is quite expected because it is energetically favourable to achieve charge neutrality 
among particles carrying conserved baryon numbers |1J . No such conservation law is followed 
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by leptons or mesons. But as soon as K~ condensation sets in, lepton fractions as well as 
H~ fraction drop. It indicates that the EoS is now mainly softened by the presence of A 
hyperons and K~ condensate. This has an interesting implication on the threshold density 
of K° condensation. It is evident from Table IV, K° condensation occurs at same density 
point 6.16 no with and without the inclusion of in our calculation. Also, a symmetric 
matter of n^,^ - and K° condensate emerges here after the onset of K° condensation. At 
a much higher density 6.7no, S° appears in the system. 

Pressure (P) is plotted against energy density (e) for various compositions of neutron 
star matter in Fig. 3. Here equations of state (EoS) are calculated with Groningen set. 
The dotted line stands for nucleons-only matter, while the dash-dotted line contains A 
hyperons in addition to nucleons. The presence of an additional degree of freedom softens 
the EoS appreciably. The solid lines correspond to A hyperon matter including K~ and 
K° condensate for antikaon optical potentials Ug(rio) = -120 to -180 MeV. The kinks on 
the equations of state mark the onsets of K condensation. Already, we have noted that 
K~ and K° condensation are second order phase transitions for all values of Ug(no) in 
our calculation. From Fig. 3, we find pressure increases with energy density even after 
the onset of antikaon condensation. The appearance of K~ condensate makes equations of 
state softer in all cases. The kinks at higher densities correspond to K° condensation which 
further softens the EoS. Also, the softness in the EoS is very sensitive to antikaon optical 
potential depth. Stronger the attractive antikaon interaction, softer is the corresponding 
EoS. 

The results of static structures of spherically symmetric neutron stars calculated using 
Tolman-Oppenheimer-Volkoff (TOV) equations JIJ and the above mentioned equations of 
state are now presented here. We have used the results of Baym, Pethick and Sutherland 
p?| to describe the crust of a compact star composed of leptons and nuclei for the low density 
(riB < 0.001/m -3 ) EoS. In the mid density regime (0.001 <n B < 0.08/m- 3 ) the results of 
Negele and Vautherin |58] are taken into account. Above this density, an EoS calculated in 
DDRH model has been adopted. The maximum neutron star masses (M max /M Q ) and their 
central densities [u cent = n cent /n ) for various compositions of matter are listed in Table 
IV. The values recorded within the parentheses correspond to the calculations including H 
hyperons in addition to A hyperons. The maximum mass of nucleons-only star is 2.313M . 
The inclusion of A hyperons softens the EoS lowering this value to 1.7O8M . Because of 
further softening due to the inclusion of H hyperons the maximum mass is reduced to a value 
of 1.620M Q . The static neutron star sequences representing the stellar masses M/M Q and 
the corresponding central energy densities (e) are exhibited in Fig. 4 for n,p,A and lepton 
matter with K~ and K° condensate and different values of U^{tiq). The softening in the 
EoS due to the presence of K condensates, leads to further lowering in the limiting masses of 
neutron stars that too attain at much earlier central densities as it is evident from Table IV. 
The maximum mass of the star varies from 1.697M Q (for Ux(no) = -120 MeV) to 1.497M 
(for Ug(n ) = -180 MeV) because strong attractive antikaon interaction in medium produces 
more softening in the EoS. For n,p,A, lepton and K condensate matter composition, K~ 
condensation thresholds occur well inside the maximum mass stars for all values of Ug (no). 
So the star is mainly composed of nucleons, A hyperons and K~ condensate. On the other 
hand, K° condensation along with K~ condensation might be a possibility in maximum 
mass neutron stars for \Ug{no)\ > 180 MeV. 
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We also inspect the effect of E hyperons on the compact star mass sequence. Already we 
have discussed that the appearance of H~ hyperons prevents the onset of K condensation for 
Uk(71o) = -120 MeV. From Table IV, we find no K condensation occurs inside the limiting 
mass neutron stars for \U^{no)\ < 140 MeV. For these values of U^(n ), the maximum 
star mass is the same as that of the case without any antikaon condensate. On the other 
hand, K~ condensate is formed inside maximum mass stars for |[/^-(no)| > 160 MeV, but 
K° condensation in neutron stars is ruled out for all values of £7^- (no) except for antikaon 
potential depth of -180 MeV. For \Ux(no) \ > 160 MeV, we observe there is hardly any change 
in the maximum masses of neutron stars compared with the cases excluding 5 hyperons. 
Already we have noted in the discussion of Fig. 2 that the density of H~ hyperon diminishes 
with the appearance of K~ condensate for \Ux(n )\ > 160 MeV. Now we show the equations 
of state for neutron star matter with and without S hyperons for [7^- (no) = -160 MeV in 
Fig. 5. The solid and dashed line represent neutron star matter with and without H 
hyperons respectively. The EoS becomes softer in the presence of S~ hyperons, but there is 
no difference between the equations of state just after the onset of K~ condensation. This 
feature is reflected in the maximum masses of neutron stars as is evident from Table IV. 

In Fig. 6, we draw the mass-radius relationship for n,p,A, lepton matter with and 
without K condensate in DDRH model using Groningen set and different antikaon optical 
potential depths and compare it with our previous result for hyperonic matter including K~ 
condensate calculated in the RMF model |J using GM1 model and U^(n ) = -160 MeV. 
The filled circles correspond to the maximum masses of compact stars. In case of no K 
condensate, the maximum mass star has a radius 11.54 km. For Ux{no) = -120 MeV, the 
maximum mass star has a radius 11.76 km, whereas it is 11.39 km for [/^-(no) = -180 MeV in 
DDRH model. The smaller radius in the latter case may be attributed to more softening in 
the EoS due to strong attractive antikaon potential. The curve corresponding to the RMF 
calculation |J has the smallest radius of 10.9 km among all the cases considered here. 

We also investigate K~ condensation in nucleons-only and n,p,A matter using Bonn A 
set as given by Table I and antikaon optical potential depth at normal nuclear matter density 
Ur^o) — -160 MeV. For Bonn A set, a-nucleon and w-nucleon couplings are density depen- 
dent whereas p-nucleon coupling is a constant one. Here we study the EoS and structure of 
neutron stars. Unlike the situation with Groningen set, antikaon condensation in this case is 
a first order phase transition which is governed by Gibbs phase rules and global conservation 
laws as given by Eqs. (35)- (39). The EoS for nucleons-only matter with and without K~ 
condensate are denoted by dotted and solid line and those of n,p,A matter are shown by 
solid and dashed line in Fig. 7 respectively. For nucleons-only matter, antikaon condensa- 
tion occurs at energy density 397.6 MeV fm~ 3 . And the phase transition is over at 570.3 
MeV fm~~ 3 . Those two points give the extent of the mixed phase. We have pure hadronic 
phase below the lower boundary and antikaon condensed phase above the upper boundary. 
On the other hand, the lower and upper boundary of the first order phase transition to K~ 
condensate in n,p,A matter for R a \ = 0.49 and R^^ = 0.4911 (case I) are shifted to higher 
energy density 415.0 and 613.2 MeV fm~ 3 respectively, because of the early appearance of 
A hyperons. We find that the effective nucleon mass becomes negative in n,p,A matter with 
and without K~ condensate respectively at 6.60no and 6.55no, where no = 0.159/m~ 3 . We 
also perform calculation for n,p,A matter with and without K~ condensate using R a \ = 0.49 
and R w \ = 0.553 (case II). In this case, K~ condensation begins at 397.6 MeV fm~ 3 and 
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the phase transition ends at 588.4 MeV fm -3 . Here A hyperons appear in the mixed phase. 
We note the EoS with and without K~ condensate for case II are stiffer compared with 
those of case I. This may be attributed to the stronger repulsion due to the larger value of 
-R^a in case II. For case II, we also get negative effective nucleon mass in n,p,A matter with 
and without K~ condensate at 6.72 uq and 6.60 no respectively. It follows from the struc- 
ture calculation using TOV equations, the maximum masses of nucleons-only stars with and 
without K~ condensate for Bonn A set are 2.55M and 2.32M Q having central densities 5.64 
and 6.84 77, , respectively. In this case, we find the radii for neutron stars with and without 
K~ condensate are 10.88 km and 9.91 km, respectively. These values of maximum masses 
and radii in Bonn A set are smaller than those of Groningen set. For n,p,A matter with and 
without K~ condensate in both case I and case II we find that the effective nucleon mass 
becomes negative before the maximum masses are reached. This feature was earlier found by 
others for n,p,A matter without antikaon condensate p6| . Because of the behaviour of the 



parameterization of the couplings in the high density regime for Bonn A set, the repulsion 
due to uj field becomes larger than the attraction of a field. Consequently, the EoS in Bonn 
A set is stiffer than that of Groningen set. A close inspection of the parameterization of 



couplings in Bonn A set has been already suggested in Ref. [31 



V. SUMMARY AND CONCLUSIONS 

We have studied K~ and K° condensation in /3-equilibrated hyperonic matter within 
a density dependent hadron field theoretical model. In this model, baryon-baryon and 
(anti)kaon-baryon interactions are mediated by the exchange of a, u, p and 5 meson. Density 
dependent meson-baryon coupling constants are obtained from microscopic Dirac Brueckner 
calculations using Groningen and Bonn A nucleon-nucleon potential. On the other hand, 
we have considered constant meson-(anti)kaon couplings in this calculation. 

For Groningen set and the values of antikaon optical potential Ux(no) = -120 to -180 
MeV, both K~ and K° condensation are found to be second order phase transition. The early 
appearance of A hyperons delays K condensation to higher density for all values of U^(n ) 
considered here. With further inclusion of H hyperons, K~ as well as K° condensation do not 
occur at all for \Ux\(no) < 120 MeV, whereas K condensate appears after being delayed by 
5 hyperons for |^|(no) > 140 MeV. It is interesting to note that as soon as K~ condensate 
appears in the system, the density of H~ drops. It is found that antikaon condensation is 
not only sensitive to the equation of state but also to antikaon optical potential depth. 

The equations of state for different neutron star matter compositions including K conden- 
sate have been studied in DDRH model. The appearance of antikaon condensation makes 
the corresponding EoS softer. This softening leads to the reduction in maximum masses 
of neutron stars for different cases considered here. For different compositions of neutron 
star matter, it is observed that K~ condensation may occur in maximum mass stars but 
the appearance of K° is ruled out except for \U^(n )\ > 180 MeV. The neutron star with 
smallest maximum mass and radius is obtained for L^j(no) = -180 MeV in DDRH model 
with Groningen set. We also studied the structure of neutron stars for nucleons-only matter 
with and without K~ condensate in DDRH model using Bonn A set. In this case, the EoS 
including K~ condensate results in a neutron star having radius < 10 km. 

We have compared the results of DDRH model with those of RMF model with GM1 set 
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||. The qualitative agreement between the results of these two models is good. Earlier, 
it was argued that many body correlations may prevent antikaon condensation to occur 
in neutron stars |TI| , |P2fl . On the contrary, the study of antikaon condensation in DDRH 
model with density dependent meson-baryon couplings which take into account many body 
correlations, shows that it is a possibility in neutron stars. In this calculation, we have 
treated meson- (anti)kaon couplings as constant. In principle, one may consider density 
dependent meson-(anti)kaon couplings. This will introduce additional rearrangement term 
in the antikaon sector. It will be reported in a future publication. 
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TABLES 



TABLE I. Density dependent meson-nucleon couplings at saturation density are obtained from 
DB calculations using Groningen nucleon-nucleon potential in Ref. [34]. Infinite nuclear matter 
properties calculated with momentum corrected meson-nucleon vertices are binding energy E/A 
= - 15.6 MeV, saturation density uq = 0.18/m~ 3 , asymmetry energy coefficient CL aS y — 26.1 MeV, 
incompressibility K = 282 MeV and effective nucleon mass m* N /rriN = 0.554. Masses for nucleons 
and mesons are m N = 939 MeV, m a = 550 MeV, m u = 783 MeV and m p = 770 MeV. The 
parameterization of density dependent a and w-nucleon couplings for Bonn A potential is taken 
from Ref. [30,31]. The nuclear matter properties in Bonn A potential are E/A = - 15.75 MeV, 
no = 0.159/m -3 , a asy = 34.3 MeV, K = 151.3 MeV and m* N /mN = 0.642. All hadronic masses for 
Bonn A case are same as in Groningen case. The p meson-nucleon coupling is density independent 
and no 5 meson is present in Bonn A case. All parameters are dimensionless. 





9aN 


9luN 


9pN 


35N 


Groningen 


9.9323 


12.1872 


5.6200 


7.6276 


Bonn A 


9.5105 


11.5401 


8.0758 





TABLE II. Scaling factor for a and uj meson-hyperon vertices for Groningen and Bonn A 
nucleon-nucleon potential. 





RwA 


R<tA 


RcrS 


RcjH 


Groningen 


0.5218 


0.49 


0.3104 


1/3 


Bonn A 


0.4911 


0.49 


0.3343 


1/3 



TABLE III. The coupling constants for antikaons (K) to a-meson, g a K, for various values of 
K optical potential depths U^uq) (in MeV) at saturation density. The results are for Groningen 
and Bonn A nucleon-nucleon potential. 





-120 


-140 


-160 


-180 


Groningen 


0.1993 


0.6738 


1.1483 


1.6228 


Bonn A 


1.1121 


1.6609 


2.2097 


2.7585 
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TABLE IV. The maximum masses M max and their corresponding central densities 
Ucent = n C ent/no for nucleon-only (np) star matter and for stars with further inclusion of hyper- 
ons (npA (£)) are given below. The results are for Groningen set. The threshold densities for K~ 
and K° condensation, u cr (K~) and u cr (K°) where u = ns/n^ and also M max and u cent for neutron 
star matter including A hyperons at different values of antikaon optical potential depth U^(no) 
(in MeV) at saturation density are given. The values in the parentheses are those of neutron star 
matter including S. 





U R (n ) 


u cr (K ) 


u cr (K°) 


Ucent 


M max /M Q 


np 








5.11 


2.313 


npA(S) 








5.13 (4.89) 


1.708 (1.620) 




-120 


3.83 (-) 


-(-) 


4.84 (4.89) 


1.697 (1.620) 


nptfA(E) 


-140 


3.17 (5.74) 


7.27 (7.39) 


4.56 (4.89) 


1.665 (1.620) 




-160 


2.65 (3.20) 


6.16 (6.16) 


4.38 (4.49) 


1.602 (1.599) 




-180 


2.28 (2.29) 


5.16 (5.16) 


5.16 (5.16) 


1.497 (1.497) 
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Figure Captions 

FIG. 1. Number densities (rij) of various particles in /3-equilibrated n,p,A and lepton 
matter including K~ and K° condensate for Groningen set and antikaon optical potential 
depth at normal nuclear matter density U^(n ) = —160 MeV as a function of normalised 
baryon density. 

FIG. 2. Number densities (n,) of various particles in /5-equilibrated n,p,A,S and lepton 
matter including K~ and K° condensate for Groningen set and antikaon optical potential 
depth at normal nuclear matter density Ux(n ) = —160 MeV as a function of normalised 
baryon density. 

FIG. 3. The equation of state, pressure P vs. energy density e for Groningen set is shown 
here. The results are for n,p and lepton matter (dotted line), n,p,A and lepton matter 
(dash-dotted line), and n,p,A and lepton matter including K~ and K° condensate (solid 
lines) calculated with antikaon optical potential depth at normal nuclear matter density of 
U R (n ) = -120, -140, -160, -180 MeV. 

FIG. 4. The compact star mass sequences are plotted with central energy density for 
Groningen set and antikaon optical potential depth of Ux(n ) = —120, —140, —160, —180 
MeV. The star masses of n,p,A and lepton matter with K~ and K° condensate are shown 
here. 

FIG. 5. The equation of state for n,p,A,S, lepton and K matter (solid line) and n,p,A, 
lepton and K matter (dashed lines) calculated with Groningen set and antikaon optical 
potential depth at normal nuclear matter density of Ux(no) = —160 MeV are shown. 

FIG. 6. The mass-radius relationship for compact star sequences for n,p, A and lepton 
matter with K~ and K° condensate for Groningen set and antikaon optical potential depth 
of U^(n ) = —120, —140, —160, —180 MeV. The mass-radius relationship for compact star 
sequence for hyperon matter including K~ condensate in the RMF model calculation (Ref. 
[9]) is also shown here. 

FIG. 7. The equation of state for n,p and lepton matter (solid line) and n,p, lepton 
and K~ matter (dotted line) calculated with Groningen set and antikaon optical potential 
depth at normal nuclear matter density of U^irio) = —160 MeV are shown. The equation 
of state for n,p,A matter with and without K~ condensate for different values of -R^a are 
also plotted here. 
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e (MeV fm" 3 ) 



